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INTRODUCTION 



New integrability case for the Riccati equation 

M. K. Malfl and T. HarkJll 

Department of Physics and Center for Theoretical and Computational Physics, 
The University of Hong Kong, Pok Fu Lam Road, Hong Kong, P. R. China 

A new integrability condition of the Riccati equation dy/dx — a{x) + b{x)y + c{x)y^ is presented. 
By introducing an auxiliary equation depending on a generating function f{x), the general solution 
of the Riccati equation can be obtained if the coefficients a{x), h{x), c{x), and the function f{x) 
satisfy a particular constraint. The validity and reliability of the method are tested by obtaining the 
general solutions of some Riccati type differential equations. Some applications of the integrability 
conditions for the case of the damped harmonic oscillator with time dependent frequency, and for 
■ " , solitonic wave, are briefly discussed. 
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\ The Riccati equations, of the type 
, ^ =aix) + bix)y + c{x)y^, (1) 

i' where a, 6, c are arbitrary real functions of x, with a,h,c G C°°{I), defined on a r eal interval / C 3? (T], Q, find 

1 . surprisingly many applications in physics and mathematics. For example, supersymmetric quantum mechanics 0], 
' variational calculus [j] , nonlinear physics [sj, renormalization group equations for running coupling constants in 

quantum field theories @, thermodynamics Q, the formulation of Newton's law the dynamical rate equations 
^ ' for physical processes driven by a combination of diffusive growth [To| . obtaining analytical solutions to the (3+1)- 
\ dimensional Gross-Pitaevskii equation in the presence of chirp and for different diffraction and potential functions [ll| , 
^ ■ and the study of the cubic nonlinear Ginzburg-Landau equations [l2| are just a few topics where Riccati equations 
^ ' play a key role. One of the main reason for their importance in physical applications is that a change of function turns 
\^ . the Riccati equation into a linear second-order differential equation, that stands as basic mathematical background for 
' many areas of physics. Group theoretical methods, which are very useful for a better understanding of the properties 
. of the Riccati equation, and a discussion of the integrability conditions from a group theoretical perspective can be 
found in [l^l . From a group theoretical point of view the nonlinear superposition principle also arises in a simple way. 
] Since the Riccati equation is a widely studied nonlinear equation, knowing that the physical system under consider- 
■ ation can be brought into Riccati form has certainly many advantages in the investigation of its properties. If one or 
' two particular solutions yf(x), i = 1, 2 of the Riccati equation are known, the equation may be solved by quadratures 

' m- 

A number of solutions of the Riccati equation can be obtained by assuming that the coefficients a{x), b(x), and c{x) 
satisfy some particular constraints. Thus, if the coefficients a, b, c, defined and continuous in some interval / C 3?, 
are related as 



H 

' , d ^ a a - P . ^ . 

a + b + c^ —log- — {ac - /3a) , 2 

ax p ap 

with a{x) and I3{x) properly chosen functions diffcrentiable in /, such that a/3 > 0, then the Riccati equation is 
integrable by quadratures [IJ]. If c{x) = 1 and the functions a{x) and b{x) are polynomials satisfying the condition 

db(x) 

A = b'^{x)- 2^-^ -4:a(x) = constant, (3) 
ax 

then 



b{x) - 

yii^)^-- ^ (4) 
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and 

6(x) - 

2/2 (.t) = (5) 

arc both solutions of the Riccati equation [l|, 0. Note that the choices c(x) = 1 and h{x) and c{x) being of polynomial 
form arc restrictive conditions of the method [ll, 01- If we know three particular solutions j/f (x), i = 1, 2, 3, then the 
Riccati equation can be solved without quadratures [11, [2[ . 

It is the purpose of the present paper to introduce a new integrability case for the Riccati equation. By introducing 
an auxiliary equation depending on a generating function /(x), the general solution of the Riccati equation can be 
obtained if the coefficients a(x), h{x), c{x), and the function f{x) satisfy a particular constraint. The validity and 
reliability of the method are tested by obtaining the general solutions of some Riccati type differential equations. 
Some applications of the integrability conditions for the case of the damped harmonic oscillator with time dependent 
frequency, and for solitonic wave, are briefly discussed. 

The present paper is organized as follows. The integrability condition for the Riccati equation is obtained in 
Section|lIl Some specific examples of the integrability condition obtained by fixing the functions b(x), c{x) and f(x) 
are shown in Section Hill Some integrable Riccati equations with fixed a(x), h{x) and f{x) are presented in Section HVl 
The physical applications of the method arc briefly outlined in Section |Vl We conclude our results in Section IVII 



II. THE INTEGRABILITY CONDITION FOR THE RICCATI EQUATION 



From an algebraic point of view Eq. ([T]) is a quadratic equation in y. We consider that its particular solutions y^{x) 
can be represented as 



-6(x) ± 



2c {x) 



where we have introduced the new function f{x) £ C°°{I), defined as 



fix)=b^ ix)-Acix) 



a {x) 



dy_ 

dx 



(6) 



(7) 



The requirement that y^{x) as defined in Eq. ([S]) is a solution of the Riccati Eq. ([T]), restricts the expression of a{x) 
to the form 



a{x) 



d 

dx 



-b{x) ± 

2c{x) 



b\x)~f{x) 
4c(x) 



(8) 



By substituting a{x) given by Eq. dH) into Eq. ([J), we obtain an auxiliary Riccati equation of the form 

h\x)-f{x) 



dy d 
dx dx 



2c(a;) 



Ac{x) 



b{x)y + c{x)y , 



(9) 



where f{x) is a solution generating function to the auxiliary Riccati Eq. Q. Therefore we obtain the following 
Theorem The general solution of Eq. ([9|) with particular solutions given by Eq. ([6]) is represented by 



2c(a;) 



(10) 



where C± is an arbitrary integration constant. 

Hence by fixing the functional forms of the functions b{x), c{x) and f{x), we can obtain the general solution of the 
Riccati Eq. ^ with a{x) given by Eq. ([8]). Alternatively, by fixing the functional forms of the functions a(x), b{x) 
and f{x), we can obtain the general solution of the Riccati Eq. with c{x) given by Eq. ([S]). 

In the case c{x) = 1 and f{x) = A = constant, from Eq. ^ we reobtain the integrability condition given by Eq. ([3]). 



III. GENERATING SOLUTIONS OF THE RICCATI EQUATION BY FIXING b{x), c{x) AND /(x) 

A. Solutions with f{x) = 

In the case f{x) = Eq. ^ becomes 



a{x) 



dx 



2c{x) 



b'^jx) 

4:C{X) ■ 



(11) 



and from Eq. (fTO|) it follows that the Riccati Eq. ([9]) satisfying the condition given by Eq. has the general solution 



- / c{x)dx + C 



-b{x) 



2c{x) 



(12) 



Thus by giving the functional form of the two functions b{x) and c{x) the general solutions of the Riccati equation 
can be obtained. 



1. Example 1. 

The coefficients in the Riccati equation 

ax A ^ 
where a, /?, ri, m are arbitrary real constants, satisfy Eq. (jlip . The equation has the particular solution 



yP{x) = -l.^KP-o^)^rf^-^ 



1 

and the general solution of Eq. ([T3|) . which follows from Eq. ([T2l) . is given by 

^^'^^ ~ a;"+i (-ax-)"^"" r(l + n, ^ax) + G ^ 2^ 
where r(a, z) — t°-~^e~*dt is the incomplete gamma function [Tsj . 



(13) 



(14) 



(15) 



Example 2. 

The Riccati equation with Bessel Jn{x) function coefhcients, given by 

dy _ ^-"+'3-1 {2xJm-iix)Jn{x) + Jrnjx) [2(-?7i + n + a - /j) J„(a:) + x (x^ J„(a;)^ - 2J„_i(x))] } 

xr 

x^ J„ (x) y + x" Jm (x)?/^ , 
where a, /3, to, and n, are arbitrary real constants, has the particular solution 

x^-"J„(x) 



yP{x) 



2J™(x) 



The general solution of Ec^. ([T6| is given by 



x^-"J„(x) ^ 1 

2Jm(a:^) C- 2-'"-ix'"+"+ir [i(m + a + l)] 1F2 [^(to + a + 1); m + 1, i(?7i + a + 3); 



(16) 



(17) 



, (18) 



where T{z) = *dt is the Euler gamma function, and pFq{a;b;z) = pFq {a;b; z) / [T (bi) ....T {bq)] is the 

regularized generalized hypergeometric function [isj . 
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B. Solutions with f{x) = constant 
By choosing f{x) = = constant, with K G 3?, Eq. (|8]) becomes 



a{x) 



dx 



-b{x) ± K 



2c{x) 



h\x) - 
Ac{x) ' 



(19) 



and from Eq. (fTO)) it follows that the general solution of the Riccati Eq. ^ satisfying the condition given by Eq. ([19]) 

is represented by 



y±{x) = e 



±Kx 



-1 

+ 


'-b{x)±K' 




2c{x) 



(20) 



1. Example 3. 



The Riccati equation 
dy 1 



dx 4 



e "^a;-"-! <i e^'^a;" [-2m + 2n + x (e'^^a;™ + 2a - 2/3)] - T 2X(n + ax) ) + 



has the particular solution 



y±{x) ^ ^e~°"'x-" {±K - e^^x") 



The general solution of Eq. (PT|) is given by 



1 



y±(x) = -a;-"e-"^ (±K - e^^z") + 



[-x(±ii: + a)]"""' r [n + 1, -a;(±A' + a)] + C± 



(21) 



(22) 



(23) 



2. Example 4- 



The Riccati equation 



^ — -e csc(na;) < e'''^ 
dx 4 



sin(TOa;) ^2q; — 2/3 + 2ncot(na;) + e'^^ sin(TOa;)J — 2TOC0s(TOa;) 
K [±K + 2a + 2ncot{nx)] > + e^"^ sm{mx)y + e"^ sin(nx)y^, 



(24) 



has the particular solution 



y±{x) = -e csc{nx) [±K - e^"" sin(r 



The general solution of Eq. (p4|) is given by 



y± 



{x) = csc(rtx) [±K - e'^^ sin(ma;)] + 



C+ - 



ef='=^+")^[(±/4r+a) sin(na;) — n cos(n2;)] 



(25) 



(26) 



C. Solutions with arbitrary function /(a 



We choose f{x) — that is, b{x) — ±^/f{x). In this case the Riccati Eq. ^ becomes 

dy 



dx 



b{x)y + c{x)y 



with the general solution obtained with the use of Eq. (fTO)) given by 



y{x) 



^Jb{x)dx 



The same result can be obtained by directly solving the Bernoulli type Eq. p7)) . 



1. Example 5. 



The Riccati equation 

dy ^ i^-n-l 

dx 4 

has the particular solution 



[-2(771 - n)x"' + a;2"+i - + (1 - 2n)^] + x'^y + x"y% 



corresponding to f{x) = x. The general solution of Eq. (|29p is given by 

2x3/2 



1 



= + 77771 



e 3 



2(n+l) 2a3/2 
3 ' 3 



C 



^. Example 6. 



The Riccati equation, 



4' 



- + 2m - 2n) + x2™+^ + (m - n)x('"+")/2 



+ .T™y + x"y^. 



has the particular solution 



yP{x) = 



^(m+Ti)/2 _ 



corresponding to /(x) = 6(x)c(.t) = .t™+". The general solution of Eq. ([5^ is given by 



2/(2^) 



1 



^(in+n)/2 _ 



C - } e "+"+2 dx 



IV. GENERATING SOLUTIONS OF THE RICCATI EQUATION BY FIXING a{x), b{x) 
Integrating Eq. ^ we obtain 

1 



c±ix) = -/±(x) [b{x)T^/m 



— / a{x)I±{x)dx + k± 
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where k± is an arbitrary integration constant, and 



I±{x) — exp 



dx 



With this form of c(x) the Riccati Eq. ([T]) can be written as 



dy 



= aix) + b{x)y + -I±ix) \ b{x) T V/M 



a{x)I± {x)dx + k± 



Then the general sohition to Eq. ([57)) is 



(36) 



(37) 



-i//±(x) b{x)T^7{x) [- j a{x)I±{x)dx + k±\ 'e±/V^'^"dx + C± 



— I a{x)I±{x)dx + 



(38) 



3. Example 7. 



The Riccati equation 



I3^x- 



rl^ o-m y 13 



—y 



(39) 



where /s e 3? is a constant, whose coefficients satisfy the condition given by Eq. ([35]) with f{x) = 0, has the particular 
solution 



2a fix^-™ 
yP(x) = -— - fce2(i-™) . 



The general solution of Eq. ([M]) is given by 



2a ^xi-"' 



4(to — l)ax'^' 



2(m- 1) 



2Ca + /31og 2a + fc/^e^u^ 



(40) 



(41) 



a;™ + (3^x 



The Riccati equation 



4. Example 8. 



da; a; 2 



(42) 



whose coefficients satisfy the condition given by Eq. psp with f{x) = 1, a(x) = /_|_^(x), and k = m, respectively, has 
the particular solution 



n/ \ mil x/l 

y (x) — —mx ' e ' . 

The general solution of Eq. is given by 

y(a;) = x"/^e^/2(x-m) + - 



^x-"^/^E^/2+i{~x/2) + C' 
where En[z) is the exponential integral function En{z) = e^^'^/t^dt. 



(43) 



(44) 



7 



5. Example 9. 



The Riccati equation 



dy (3 a 



y + 



2j3x +{a + l)kx- 



dx X X' 2 

whose coefRcients satisfy the condition given by Eq. ([55)1 with f{x) ~ x"'^ , has the particular solution 

[a + \)kx^ +2(3 



yP{x) 



1 - a2 



The general solution of Eq. p5|) is given by 

y{x) = — - kx 2 

a + 1 



A(3x 



(a2 _ l)x2i^l 



Q+3 . (a+l)kx 2 



a+l ' ' a+1 ' 



2/3 



4C/3 



where 2^1 (a, b] c; z) is the hypergeometric function. 



(45) 



(46) 



(47) 



V. APPLICATIONS IN PHYSICS 

In the following we consider two physical application of the new integrability conditions for the Riccati equation 



A. The damped time-dependent harmonic oscillator 



Many natural processes can be modeled in the classical regime by the motion of a damped particle in a time- 
dependent harmonic potential, described by the equation [l6j 



x + "f{t)x + uj'^{t)x = 0, 



(48) 



where a dot denotes the derivative with respect to the time t, x{t) is the position of the particle, ^(t) is the damping 
function, and a;^(i) is the time dependent angular frequency, respectively. Riccati parameter families of damping 
modes, related to the Newtonian free damping ones by means of Wittens supersymmetric scheme were considered in 
[isj . This procedure leads to one parameter families of transient modes for each of the three types of free damping, 
corresponding to a particular type of anti-restoring acceleration. The Ermakov-Lewis procedure was applied to the 
one-parameter damped modes introduced in [l^ in [l^ . 

By introducing the transformation u ~ i/x, x{t) = xq cxp j^u{t')dt' , where xq = x(0), Eq. (j48p can be 

transformed into a Riccati type equation, given by 



(49) 



Eq. (Uni) has three general solutions, corresponding to three different constraints imposed on the damping and 
angular frequency functions. If the two functions uj(t) and "f{t) satisfy the condition 



Mt) , i\t) 



the general solution of Eq. P^]) is given by 



7 ^ 1 

uit) = — ^ H . 

^ ' 2 C + t 



(50) 



(51) 



The integration constant C can be determined from the initial condtion u(0) = x(0)/x{0) — Vq/xq^ where Vq is the 
initial velocity of the particle. Therefore the integration constant is 



C 



2xq 



2vo + 70x0 



(52) 
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where 70 = 7(0). The general solution of Eq. (|48|) is given by 



x{t) = xo ( 1 + — ) cxp 



7 (f) dt' 



If the condition 



2 4 

is satisfied, where K is & constant, the general solution of Eq. (|49p is given by 



u±{t) = - 



7(t) ± K 



exp {±Kt) 



C± ± {l/K)cxp{±Kt)' 



and the general solution of Eq. (|^5)) is given by 



x± (t) 



Xo 



(C± ± JT-i) 

For this case the value of the integration constant C± is 

2x±a 



_e-(i/2) /o'' -/it')dt' ^ /^-ie±-ff*/2^ _ 



C± = 



2v±o + (70 ± K) x±o 
where u±o = v±o/x±o. Finally, if there is a function f(t) so that the condition 

iHt)-fit) 



Jit) ± 



+ 



holds for all t, then the general solution of Eqs. (|49|) and (|48| are given by 



7W ± 



and 



a;-|-(t) = x±Q exp 



u± (t') dt' 



(53) 



(54) 



(55) 



(56) 



(57) 



(58) 



(59) 



(60) 



respectively, where x±o is an integration constant. 

In order to obtain the values of the arbitrary integration constants for the solutions given by Eqs. ([5^ and ([50)1 
from the initial conditions a;±(0) ~ x±o and i±(0) = w±0: smd from the initial values of the functions f{t) and j{t), 
we denote / ^Jf{pjdt = F{t), and j e^^^^^dt — ±G{t), respectively. The initial values of these functions a.t t — 
are F{0) = Fq, and G{0) = Gq. The t = value of /(t) is /(O) = /o. Then the values of the arbitrary integration 
constants C± can be obtained as 



C± 



a;±o 2 



TGo. 



(61) 



B. Solitons 



Second order partial differential equations of the form 



b [x, t) ^-j^ + a (x, t) ^^A^ + V [x, t) * [x, t) = 0, (62) 

OX'^ VOL 

where z; is a constant, and a, 6, V are arbitrary functions of the coordinates x and are used for the description of a 
large variety of physical models. Of particular importance are the so-called soliton (solitary wave) solutions pol [2l|. 
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in which the dependence of the wave function ^ {x,t) is assumed to be of the form ^' {x,t) = {vt — x) = ^ (^). 
Therefore sohtons can be obtained as solutions of the second order ordinary differential equation 

b (0 *" (0 + « (0 iO + viO^ iO = 0, (63) 

where a prime denotes the derivative with respect to ^. By dividing this equation with b{^) ^ 0, and denoting 
= a(^) /b{C) and (^) = V /b{^), Eq. §3^ takes the form of Eq. dH]), the equation of motion of the 
damped time dependent harmonic oscillator, 



(0 + 7 (0 (0 + (0 * (0 = 0. (64) 

Therefore all the results of the previous subsection can be applied to the case of the solitons, leading to the possibility 
of constructing explicit exact solitonic solutions for the wave-type equations of mathematical physics. 



VI. CONCLUSIONS 



A new method to generate analytical solutions of the Ricacti equation was presented. The method is based on 
the correspondence between the initial Riccati equation and a more general equation containing a solution generating 
function f{x). If the coefficients of the Riccati equation and the function f{x) satisfy one differential and one integral 
constraints, the general solution of the Riccati equation can be explicitly obtained. The method was illustrated by 
obtaining the general solution for a number of specific Riccati type equations, with coefficients satisfying the required 
integrability condition. Some physical applications of the integrability method were considered, by explicitly obtaining 
several classes of general solutions for the harmonic damped oscillator with time dependent frequency. The possibility 
of obtaining soliton type solutions of the second order partial differential equations was also briefly considered. 
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